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The third interaction virial coefficient Caww for air-water vapor mixtures is estimated in the tem- 
perature range to 100 °C, by means of an approximate method based on molecular association. The 
results are believed accurate to within a factor of two. The value of Caww can be combined with previ- 
ous estimates of the other interaction coefficient Caaw and the values of the third virial coefficients of 
pure air and pure water vapor, to give a quantitative description of the equation of state of air-water 
vapor mixtures through the complete third virial coefficient. It is shown that the maximum total con- 
tribution from both Caww and C a <iw to the compressibility factor is only of the order of parts in 10 4 up to 
100 °C and 100 atmospheres total pressure. At —50 °C, it can be shown by extrapolation of Caww that 
the contribution of the C a ww term to the compressibility factor is only on the order of parts in 10 11 . 
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1 . Introduction 

The thermodynamic properties of air-water vapor 
mixtures near atmospheric pressure have been deter- 
mined and tabulated by Goff and his co-workers [1-4]. * 
Close correlation was found between the equation of 
state in the virial form and experimental observations; 
since the virial form has a sound theoretical basis, 
Goff was able to make accurate interpolations and 
extrapolations. In this work, Goff neglected terms 
involving the third interaction virial coefficients C aa w 
and Caww, which account respectively for interactions 
among two air molecules and one water molecule, and 
one air molecule and two water molecules. This 
neglect was dictated at least in part by necessity, since 
no information on C a aw and Caww was available at that 
time, but the good agreement between theory and 
experiment indicated that the neglect was probably 
justified. However, the effects of C a aw and C a ww in- 
crease with increasing pressure, and it cannot be taken 
for granted that they can be neglected at high pres- 
sures, where interactions among three molecules 
become important. The situation was improved 
somewhat in 1963, when Mason and Monchick [5] 
presented theoretical estimates of Caaw, believed accu- 
rate within a factor of two. Only Caww then remained 
unknown in the treatment of the equation of state of 
moist air through the third virial coefficient. 

The purpose of this paper is to present theoretical 
estimates of Caww- Internal evidence in the calcula- 



tions indicates that the results are reliable within a 
factor of two from to 100 °C; even order of magnitude 
reliability would be valuable, however, because no 
other information whatever is available on Caww The 
results can be extrapolated down to —50 °C with no 
loss in overall accuracy in the equation of state, since 
the maximum water content of moist air is so small at 
low temperatures. The results can also be extrapo- 
lated above 100 °C, but the upper limit of usefulness 
is uncertain. At 100 °C and pressures up to 100 atm, 
however, the total contribution from both Caaw and 
Caww to the compressibility factor is only of the order 
of parts in 10 4 , so considerable extrapolation is possible 
if only moderately accurate equation of state informa- 
tion is needed. 

The rigorous calculation of C a ww is difficult because 
one of the pair interactions, ww, is qualitatively dif- 
ferent from the other two, aw, inasmuch as a long- 
range dipole-dipole force is involved. However, if 
the dipole-dipole interaction is strong enough, a good 
approximation can be obtained by assuming that the 
two water molecules form a quasi-stable dimer, w; 2 , 
whose interaction with a single air molecule is rela- 
tively easy to calculate. This is the physical basis for 
the present theoretical estimates. 

2. The Virial Equation of State 

The virial equation of state through the third virial 
coefficient for a mixture is 
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which can also be written as a series in the pressure, 



PV 



RT 

where 



1+fln 



RT 



+ (C mix — fiLv) \~^) 



(2) 



P = absolute pressure, atm, 
V= molar volume of mixture, cm 3 /mol, 
T= absolute temperature, °K, 
R = gas constant, cm 3 -atm/mol-°K, 
^mix = second virial coefficient of the mixture, 

cm 3 /mol, and 
C mix = third virial coefficient of the mixture, (cm 3 / 
mol)' 2 . 

The virial coefficients of a mixture depend on both 
temperature and composition, but the composition 
dependence is simple. For a binary mixture of a and 
w , it is 

** mix = % a aa ' -^ a ^ w "aw ' ^iJ^ww \o) 

^ mix = % a^aaa " <->X ^Cw^>aaw ' *->% a%> w ^ aww ' X w L>www (4) 

where 

#a, #u> — mole fractions of air and water vapor, 
respectively, 

Baa, B ww = second virial coefficients of pure air and 
pure water vapor, respectively, 

B aw = second interaction virial coefficient for 
one air molecule and one water 
molecule, 

Caaa, C www = Third virial coefficients of pure air and 
pure water vapor, respectively, and 

Caaw , Caww = third interaction virial coefficients for 
two air molecules and one water 
molecule and for one air molecule 
and two water molecules, respec- 
tively. 

All the Bij and Ctjk are functions of temperature only, 
and are independent of composition. 

3. Association Treatment of Virial Coefficients 

Our approximation is based on the treatment of the 
equation of state in terms of chemical association, as 
has been discussed, for example, by Woolley [6]. In 
this treatment, a single-component gas is considered 
to be an ideal mixture of molecular clusters consisting 
of monomers, dimers, trimers, etc. The actual non- 
ideal behavior of the gas is thus accounted for by the 
change in the total number of molecules due to the 
association — an ideal gas would consist only of 
monomers. This treatment leads to the following 



equation of state for one mole of a single-component 
gas: 



j^==l-K 2 P + (3Kl-2K 3 )P*+. 



(5) 



where K 2 is the equilibrium constant for dimer forma- 
tion, K3 that for trimer formation, and so on. 

At this point the results are completely rigorous, and 
in fact there is a one-to-one correspondence between 
the equilibrium constants and the Mayer cluster inte- 
grals of statistical-mechanical theory. Approximations 
enter only in the evaluation of the equilibrium con- 
stants, and the chemical association method is most 
useful in those cases where the interactions are strong 
enough that K 2 , for instance, can be obtained by 
considering only those pairs of molecules which are 
bound together in their lowest energy state. If all 
states of binding are considered, then the ordinary 
virial equation of state is recovered, but couched in 
somewhat different language than usual. 

For example, for pure water vapor, treated as a 
mixture of monomers, dimers, etc., comparison of 
eq (5) and the virial equation of state for water monomer 
shows that 

B WW = -K 2 RT, C www = (4Kl-2K 3 )(RT)\ etc. (6) 
According to the law of mass action, 



K 2 = 



x%RT 



(7) 



where x W2 is the mole fraction of dimers and x w the 
mole fraction of monomers, 
and (7) leads to 



Combination of eqs (6) 
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(9) 



for x w ^ x w - 

The treatment is easily extended to gas mixtures, 
and a hybrid description is even possible, in which 
strongly interacting components of the mixture are 
handled by the chemical association method, and 
weakly interacting components by the ordinary virial 
coefficient method. Such a hybrid description is 
appropriate for mixtures of air and water vapor. We 
will not, however, calculate B ww from eq (6), but will 
take B ww as an experimental quantity and use eq (8) 
to calculate the mole fraction of dimers, which will 
be needed in the estimation of C a ww 

Barua, et al. [7] have shown that for polar molecules 
a fairly good approximation to B and K 2 can be ob- 
tained by considering only those pairs of molecules 
which interact strongly enough to be bound together 
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in the lowest energy state — that is, by considering 
only dimers whose component molecules have their 
dipoles alined head-to-tail. The contributions from 
other dimer configurations can be ignored if the tem- 
perature is not too high. This study thus gives quan- 
titative evidence for the validity of the concept of 
dimers of water molecules. 

For a description only through the third virial co- 
efficient, we can consider the three-molecule inter- 
actions either as trimer formation, or else as the 
nonideal interaction between a dimer and a monomer. 
Our model of a mixture of air and water vapor through 
the third virial coefficient will be that of a mixture 
of three species — air monomers, water monomers, 
and water dimers. All of these species interact with 
one another except for two water monomers, whose 
interaction is already taken into account by the ex- 
istence of the water dimers. The equation of state 
of this ternary mixture is 

PV 1 

^f = l + « 2 #«« + 2«Baw + x'^B ww ) - 



i (Zx a X W2 Daw 2 '~^ w Xw2^ww2) 



I 



+ « 3 C aaa + 3< 2 <C a ^)^+. . . (10) 

in which we have written B ww explicitly rather than 
use eq (8), since we will take B ww as an experimental 
quantity. Other interactions, such as those involving 
B W2 w 2 or Caaw 2 , contribute to the fourth virial coefficient 
and are neglected. The mole fractions have been 
written with primes to indicate that they are not the 
same as those appearing in eqs (3) and (4). 

The degree of dimerization is on the order of parts 
per thousand of the total number of water molecules. 
We therefore make the approximation that the x' a 
and x' w of eq (10) are equal to the x a and x w of eqs (3) 
and (4). By comparison of these equations, we see 
that 

^mix~ \£x>aX>Wzt>aw2 * ZXwPC>w^ww%) * 

"T \X a L>aaa ' ^X a Xw^aaw) • 1-t-U 

Solving eq (8) for x W2 and substituting into the first 
term on the right, we find 



OXaX w L>aww~r Xyf^www Zx a X ^an^^ww ZX' u p ww 2 ^ u 



and therefore 



B u 



)— B a w 2 



Zd WW9 d u 



(12) 



(13) 



(14) 



The quantity sought, C a ww, has therefore been 
related to a product of second virial coefficients, one of 
which (B ww ) is known experimentally, and the other 
of which (BawJ remains to be calculated. In addition, 
if C W ww is known experimentally, we can use eq (14) 
as a check on the internal consistency of our above 
calculations or as a source of information on the dimer 
properties (through B WW J. 

4. Calculation of B WW2 

We can now calculate values of Bww from eq (14), 
and compare with values obtained by methods de- 
scribed in this section. This affords a check on our 
procedures to this point. 

It must be emphasized that we are not trying to 
derive values for B ww as an end in itself— rather, we 
wish to determine the ranges within which the molecu- 
lar potential parameters for the dimer-monomer inter- 
action lie. This is important, for these parameters will 
be used to obtain, in order, the molecular potential 
parameters associated with the air-dimer interaction, 
Baw 2 itself, and finally C a ww from eq (13). 

The values of B ww and C www are listed in table 1, as 
taken from the correlations of Goff [4] and Keyes [8]. 
The respective values of B ww warrant using the average 
value in all calculations. The respective C www values, 
however, should be considered individually. Goff [4| 
assigns a tolerance of ± 120 percent at 60 °C and of 
±30 percent at 90 °C to his values of C W ww, with the 
tolerances taken as twice the probable error. This 
points out the large uncertainty in this quantity, par- 
ticularly when we have used GofFs relations to extrapo- 
late his values below 60 °C. 



TABLE 1. Experimental values of the second and third virial 
coefficients of water vapor; extrapolated values in parentheses 







- B, rir 




-C 








enr'/mol 




liter*/ 


mol s 
















Goff [4] 


Keyes [8] 


Ave. 


Goff [4] 


Keyes [8| 


°c 















(1832) 


1854 


1843 


(64.3) 


93.8 


10 


1510 


1520 


1515 


(34.7) 


49.3 


20 


1260 


1266 


1263 


(19.7) 


27.6 


30 


1074 


1070 


1072 


(11.7) 


16.2 


40 


924 


917 


920 


(7.18) 


9.87 


50 


803 


794 


798 


(4.51) 


6.26 


60 


705 


695 


700 


3.01 


4.06 


70 


625 


613 


619 


2.03 


2.76 


80 


558 


546 


552 


1.39 


1.92 


90 


501 


489 


495 


0.98 


1.37 


100 


(453) 


441 


447 


(.70) 


0.99 



Two sets of B WWo values are calculated from eq (14), 
one set from GofFs C W ww, the other from Keyes' C WW w 
The results are given in table 2. 

Other estimates may be obtained from assumptions 
about the intermolecular forces. We represent the 
interaction between two molecules as the sum of a 
spherically symmetric Lennard-Jones (12-6) interac- 
tion and the interaction between point dipoles em- 
bedded in the center of each molecule. The potential 
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energy <pij{r) between molecules i and j a distance r 
apart is thus 



<Pij(r) = 4,€ij 



®~-ffl\ 



Ma 



— (2 cos Oi cos Oj — sin 0, sin Oj cos </>) (15) 

r 

where /jlj and /Xj are the dipole moments, 0; and Oj are 
the angles between the dipole axes and the line of 
centers of the molecules, </> is the azimuth angle of the 
axes, and €ij and (Tjj are adjustable constants. If either 
dipole moment is zero, eq (15) reduces to the Lennard- 
Jones (12-6) potential, for which eij is the depth of the 
potential well and (p(<Jij) = 0. Equation (15) is usually 
called the Stockmayer potential [9]. 

TABLE 2. Experimental and calculated values of the water monomer- 
dimer second virial coefficient 









-flrr.l 


ter/mol 






Temp. 
















Exptl.. eq 


(14) 


Calc. 


Calc. 


Calc. 


Calc. 




(Keyes'C u . u . u )(G 


ofTs Cuuu) 


eq (18) 


eq (19) 


eq (21a) 


eq(21b) 


°C 

















25.4 


17.4 


21.53 


11.91 


3.39 


5.08 


10 


16.3 


11.5 


17.52 


9.84 


2.66 


3.99 


20 


11.0 


7.80 


12.96 


8.20 


2.25 


3.15 


30 


7.56 


5.46 


10.44 


6.72 


1.89 


2.56 


40 


5.36 


3.90 


8.21 


6.13 


1.62 


2.20 


50 


3.92 


2.8 r 


6.61 


4.84 


1.37 


1.82 


60 


2.90 


2.k 


5.58 


4.22 


1.21 


1.55 


70 


2.23 


1.64 


4.27 


3.67 


1.09 


1.36 


80 


1.74 


1.26 


3.65 


3.16 


0.96 


1.18 


90 


1.38 


0.99 


3.10 


2.73 


.86 


1.06 


100 


1.11 


.79 


2.63 


2.46 


.79 


0.95 



The values of etj and aij must be determined from 
some experimentally measured property. For pure 
air, which has no dipole moment, the values used by 
Mason and Monchick [5] and by Chaddock [10] are 
<T aa = 3.522 A and e a Jk = 99.2 °K, where k is Boltz- 
mann's constant. For pure water vapor, Rowlinson 
[11] gives 0-^ = 2.65 A and €u^/A: = 380 °K from 
analysis of B ww , and Monchick and Mason [12] give 
<r ww = 2.1\ A and € ww /k = 506 °K from analysis of 
viscosity data. Tabulations of virial coefficients for 
the Stockmayer potential are given in terms of a 
reduced dipole parameter, defined as 



« = 



(Mjtl 



V8 



(16) 



tijdtj 



The value of t% w is 1.2 according to Rowlinson and 0.85 
according to Monchick and Mason, with the dipole 
moment of water taken as /jl w = 1.83 X 10~ 18 esu. 

Stogryn and Hirschfelder [13] have shown how to 
estimate €ij and crij for dimer-monomer interactions 
from those for monomer-monomer interactions, accord- 
ing to simple geometrical arguments. The results 
have been elaborated by Barua and Das Gupta [13] on 
the basis of the pressure dependence of the viscosity 
of superheated steam. They recommend the values 



€ww 2 — I.oZEwwi 0"ww 2 LlOCTti 



(17) 



monomers with their dipoles in a head-to-tail configura- 
tion, the dipole moment of the dimer is twice that of 
the monomer, ll w = 3.66 X 10~ 18 esu. 

' 2 

Once e, o% and t* are known, it is straightforward 
to calculate B W u 2 from published numerical tables 
[9, 11] for the Stockmayer potential. 

The Rowlinson [11] parameters for water vapor, 
based on virial coefficient data, are considered first, 
since so far only virial coefficient data have been used 
in this presentation. 

Equations (16) and (17) lead to 

t$w 2 =1.2, cr ww =3.07 A, €«*/* = 501 °K. (18) 

The values of B ww obtained from these parameters 
are found in table 2. The agreement between the 
values of —B aw calculated from eqs (14) and (18) is 
generally good, although the latter are somewhat 
higher except at °C. The implication is that the 
assumptions made to derive our model are at least as 
reasonable an approximation of the true physical 
situation as those used in the derivation of the Stogryn 
and Hirschfelder [13] model. 

If the viscosity-derived parameters of Mason and 
Monchick [5] for water vapor are used with eqs (16) 
and (17), we find 



tU =0.83, a-ww =3.14 A, 6ww Jk = 66S °K. 



(19) 



Since we assume that the dimer consists of two 



The values of B ww derived from these parameters are 
also given in table 2. The curve of this set of B WW2 
values versus T is flatter at the low temperature end 
than the three curves obtained above from virial coeffi- 
cient data, but the agreement still is not unreasonable. 

Since the values of B ww derived from eq (14) do 
agree fairly well with those derived from other sources, 
we now seek, by curve fitting, improved Stockmayer 
potential parameters which will describe the values 
from eq (14). These will be needed in the calcula- 
tion of Caww 

The curve-fitting should be done so as to produce 
preferentially better agreement at high temperatures 
than at low temperatures for two reasons. First, the 
experimental values of C www are less uncertain at high 
temperatures. Second, and more important, is the 
fact that the low-temperature values of B WW2 are dom- 
inated more by t^ Wz than are the high-temperature 
values, but the parameter t% W2 is of no interest for the 
present purposes, since it does not contribute at all to 
the determination of Caww In other words, we want 
the derived values of e WW2 and cr WW2 to represent the 
spherical (12-6) part of the potential, not the dipole 
part, as there is no dipole contribution to B a w 2 - 

Therefore, we require that the B wv>2 values obtained 
by fitting should lie between —0.79 liter/ mol and 
— 1.11 liter/mole at 100 °C, that is, between the values 
predicted respectively by using GofFs and Keyes' C www 
in eq (14). At the same time the slopes of the curves 
in the 100 °C region must be almost the same to prevent 
any gross divergences at higher temperatures. 



222 



The limiting ranges of cr W w 2 and e WW2 are set arbi- 
trarily, but not unreasonably, in the following manner. 
The lower limits are taken to be the corresponding 
values for pure water vapor given by Rowlinson [11]. 
This is equivalent to setting the coefficients of cr ww and 
€ wv in eq (17) equal to 1, a change of 13.8 percent in 
the former and 24.2 percent in the latter. The upper 
limits are set by raising the coefficients by the same 
percentages, and using the cr ww and e W w parameters 
of Mason and Monchick [5] in eq (17). The limits on 
t ww are found by using the limiting values of e W w 2 and 
a WW2 in eq (16). 

This procedure leads to 



(crww) 



ww) Rowlinson 



= 2.65 A, 



Gwwjk^ \€wwlk) Rowlinson ^oU K, 
l ww 2 ^ ^ ^* 

cr W w 2 ^ l-32(cr u , M ,) M and M = 3.58 A, 

€ww 2 lk ^ 1.64 (e ww lk) M . and m. = 830 °K, 



t* ^ 0.45. 



(20a) 



(20b) 



ancy is open to some speculation, but we remind the 
reader again that first, we demanded a fit at the high 
temperatures to have the cr and € values be more repre- 
sentative of the spherical part of the potential, and 
second, we are seeking various estimates of cr WW2 and 
€ww 2 rather than attempting to derive B WW2 . 

Summarizing this section, we have obtained four 
sets of cr ww and e WWl parameters (t* is of no further 
interest). Two were obtained through the method of 
Stogryn and Hirschfelder [13], using first experimental 
virial coefficient data and then experimental viscosity 
data, and two from curve-fitting, subject to several 
restraints, to the B ww values derived with the methods 
of this paper. We also obtained what are thought to 
be reasonable outside limits on the values of cr ww and 

€ww 2 ' 

5. Calculation of C aW w 

The various values of cr ww and € WW2 derived in the 
previous section may now be used to obtain the param- 
eters daw and e a w 2 for the air-dimer interaction, from 
which Bqw and finally C a ww may be obtained. The 
procedure is based on semi-empirical combining rules 
[9] for € and <r, 



Within these ranges, the predicted curves are slightly 
flatter in the 100 °C region than those predicted by eq 
(14). This means that a curve which predicts —1.11 
liter/mole at 100 °C will diverge immediately from 
higher temperature B ww values derived from Keyes' 
C www . In order to prevent this immediate divergence, 
we changed the lower permissible B ww „ limit at 100 °C 
ii om —1.11 liter/nwi to —0.95 liter/ mol, that is, miawav 
between the limits set by eq (14). 

Empirical curve fitting with B ww matching at both 

-0.79 liter/mol and -0.95 liter/mcJlit 100 °C led re- 
spectively to these choices of parameters: 

t* W2 =l.U, o- M , M , 2 = 3.41 A, €^£ = 380 °K (21a) 

1*^=1.25, o-^ 2 = 3.30 A, e^/A—380 °K. (21b) 

The cr ww values are not unreasonable. The fact 
that e ww Jk is the same as the value for pure water 
vapor (lower limit) in both cases is perhaps annoying, 
but can be accepted in the present context without 
undue alarm. 

The values of B* versus T* for the parameters 
were interpolated graphically between tabulated values 
of t* [9], and these in turn v^ere extrapolated graph- 
ically to the lowest temperatures. Therefore the 
derived parameters may not be the "best" ones in a 
least squares sense, but they are certainly close enough 
for our purposes. 

The values of B WW2 corresponding to eqs (21a) and 
(21b) are given in table 2, and clearly the agreement 
with the curves derived from eq (14) is poor except at 
the highest temperatures. The reason for the discrep- 



€ij = {ene^l\ try = (l/2)(<r« + <rjj). 



(22) 



These rules work fairly well for both polar and non- 
polar molecules [15]. Applying these rules twice, we 
obtain expressions for the desired parameters e a w 2 
and cr a w 2 , 

€aw 2 — €ww 2 \€aal €ww) <> 

(T(iw 2 = & W w 2 + (l/2)(o-«« — a ww ). (23) 

The values of cr aM , and e aW2 are derived from the four 
sets of wwi parameters obtained in the previous sec- 
tion. For the three associated with virial coefficient 
data [eqs (18), (21a), (21b)] we use in eq (23) Rowlinson's 
[11] values for €ww and (7ww, while for the set associated 
with viscosity data [eq (19)] we use in eq (23) the Mason 
and Monchick [5] values. We find 

€awjk = 256 °K craw= 3.51 A [from eq (18)] (24a) 
€awjk= 194 °K craw 2 = 3.85 A [from eq (21a)] (24b) 
€awjk= 194 °K ow = 3.74 A [from eq (21b)] (24c) 



zaw 2 



/k = 296 °K (Taw =3.55 A [from eq (19)] (24d) 



The last set is of particular interest, since it is 
derived entirely from viscosity data. The viscosity 
is relatively insensitive to the dipole interaction [12], 
so the values of €^2 and cr W w 2 from which eq (24d) is 
derived should be good representations of the spherical 
part of the potential. 
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Two more sets of aw>i parameters should be derived— 
namely, those associated with the limiting cases of the 
ww 2 parameters [eqs (20a) and (20b)]. From these we 
obtain respectively 

€awjk = 194 °K <Taw 2 = 3.09 A (25a) 

e aw Jk = 423 °K cr aW2 = 4.02 A. (25b) 

From each of these six sets of parameters, we cal- 
culate first B aw . z and then, using eq (13), C a ww- The 
results of the latter are given in table 3. 

Table 3. — Calculated value of the third interaction virial coefficient 
C aW w based on the parameters of eqs {24) and (25) 





- C„ ww , liter 2 /mol 2 


Temp 


Virial 


Viscosity* 


Limits 




Eq (24a) 


Eq (24b) 


Eq (24c) 


Eq (24d) 


Eq (25a) 


Eq (25b) 





0.152 


0.121 


0.111 


0.203 


0.062 


0.549 


10 


.117 


.090 


.083 


.156 


.047 


.424 


20 


.092 


.072 


.066 


.123 


.037 


.332 


30 


.073 


.057 


.052 


.098 


.029 


.266 


40 


.058 


.046 


.042 


.080 


.024 


.216 


50 


.048 


.037 


.034 


.065 


.019 


.178 


60 


.040 


.030 


.028 


.054 


.016 


.148 


70 


.034 


.025 


.023 


.046 


.013 


.123 


80 


.028 


.021 


.019 


.039 


.011 


.105 


90 


.024 


.018 


.016 


.033 


.009 


.090 


100 


.021 


.015 


.014 


.028 


.008 


.077 



*Best estimates. 

Those derived from experimental viscosity values 
[i.e., from eq (24)] are felt to be the most reliable. This 
is based on the relative independence of viscosity on 
the dipole interaction, and on the success of the dimer- 
monomer model in predicting the pressure dependence 
of the viscosity of steam [13, 14]. 

The viscosity-based values differ from the limiting 
values derived for eq (25a, b) by at most a factor of 
three, and from all other values by no more than a 
factor of two. 

We are confident that the values are accurate to 
within a factor of four, and probably to within a factor 
of two. 

6. Discussion 

We are now in a position to estimate the contribu- 
tions of Caaw and C aW w to the compressibility factor of 
moist air; that is, to estimate the magnitudes of the 
terms (SxlxwCaaw/V 2 ) and {SxaxljCaww/V 2 ), as compared 
to unity. Since the maximum water vapor content of 
air increases rapidly with increasing temperature, the 
largest contribution of these terms comes at high tem- 
peratures, even though C (l aw and C aW w decrease in 
magnitude with increasing temperature. At 100 °C 
the partial pressure of water vapor can be taken as 
1 atm, and the mole fractions as 



where P a is the partial pressure of the air. These 
relations are only approximate for nonideal gases, but 
adequate here because in C mix we are seeking only the 
highest-order deviation terms. Similarly, the molar 
volume of the mixture is 



V=RTI(P a + l). 



(27) 



Taking C aaW ; = + 696 (cm 3 /mol) 2 [5] and C aW w = — 2.8 
X 10 4 (cm 3 /mol) 2 (table 3), we find from eqs (26) and 
(27) that 

$x*x v >CaawlV 2 = 2.23 X 10- 6 Pg/(P a + 1) (28a) 

SXaXlCawwIV 2 = 0.90 X lO^PJ (P a + 1). (28b) 

Thus even for P a = 100 atm the contribution of these 
two terms is only of the order of parts in 10 4 . 

The present results are readily extendable to higher 
or lower temperatures by means of the awi parameters 
of eq (24d) in conjunction with tabulated virial coef- 
ficients for the (12-6) potential [9]. There is no loss 
in overall accuracy in PV/RT by extrapolation to low 
temperatures. For example, a calculation similar 
to that performed above, based on a temperature of 
-50 °C (P W , = 3.9X10- 5 atm) with P a = 100 atm and 
C a ww = — 1 X 10 6 (cm 3 /mol) 2 , indicates that the contri- 
bution of the Caww term to the compressibility factor 
is on the order of parts in 10 11 . There is loss in ac- 
curacy by extrapolation to high temperatures because 
of the increasing water vapor content; how much ex- 
trapolation is permissible depends on how much un- 
certainty in PV/RT is acceptable. It must also be 
remembered that the physical model of dimer forma- 
tion becomes less accurate as the temperature 
increases. 
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